Transmission phase shift of phonon-assisted tunnehng through a quantum dot 
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The influence of electron-phonon interaction on the transmission phase shift of an electron passing 
through a quantum dot is investigated by using the scattering theory. The transmission phase versus 
the intra-dot level shows a serial of phonon-induced dips. These dips are highly sensitive to electron- 
phonon interaction strength A, and they are much more pronounced than phonon-assisted sub-peaks 
appeared in the conductance. Phonon-induce dephasing is also studied, and the results show that 
the dephasing probability Td monotonically increases with the electron-phonon interaction strength 
A. The dephasing probability Td oc A^ for small A but Td oc A at large A. 

PACS numbers: 73.63.Kv, 71.38.-k, 03.65.Vf, 03.65.Nk 



I. INTRODUCTION 

Electron transport through a mcsoscopic system, e.g. 
a quantum dot (QD), has been extensively investigated 
in the last two decades. Due to the fact that the size 
of a mesoscopic device is within the phase coherent 
length, the phase of a wave function plays a key role 
on the electronic transport. So the transmission am- 
plitude t = |t|e*^, which describes the electron tunnel- 
ing through a mesoscopic system, is a complex number. 
Its magnitude square |tp is the transmission probability 
which is observable in the measurement of current or con- 
ductance. The transmission phase 9 describes the phase 
change when an electron tunnels through a device. This 
phase 9 is in general lost in the measurement of current 
or conductance, so that it is difficult to acquire in ex- 
periment. Using an AB interference ring device, Yacoby 
et alr^ tried, for the first time, to measure the transmis- 
sion phase 9 through a QD. Couple years later, Schuster 
et a& utilized an open multi-terminal AB ring device to 
successfully measure the phase 9. Since then, to investi- 
gate the transmission phase has generated a great deal of 
theoretical and experimental interest with a fare amount 
of efforts focusing in this field. On the experimental side, 
for example, Buks et ali^ reported that controlled deco- 
herence could be achieved in a device with a QD that 
is capacitively coupled to a quantum point contact in 
close vicinity. The phase evolution in the Kondo regime 
was experimentally investigated a few years back;^ and 
was found to be highly sensitive to the onset of Kondo 
correlation. Recently, Leturcq et al^ investigated the 
magnetic field symmetry and the phase rigidity of the 
nonlinear conductance in a AB ring. On the other hand, 
the success of these experiments generates a number of 
theoretical studies. In early of 1980s, Buttiker found the 
phase rigidity in a two-terminal AB ring device due to the 
time-reversal symmetry and the current conservationii 
After the experiment by Schuster et al.^ many follow-up 
theoretical efforts focused on and tried to interpret the 
measured results of the transmission phase 9 through a 
QD, in particular, the abrupt lapses of 9 between two 



adjacent resonances and the similar behavior of 9 for 
all resonant peaksi^i^ In addition, some works have also 
studied the transmission phase in the Kondo regime^ or 
with the photon-assisted tunneling process under a time- 
dependent external fieldfi^ etc. 

Another subject, the electron-phonon (e-ph) interac- 
tion in a single- molecular QD, has also generated a great 
deal of interest in recent years. The phonon-assisted tun- 
neling peaks or steps have been experimentally observed 
in various single-molecule transistor systems Ji^^iii^ Park 
et aijii observed phonon-assisted tunneling sub-steps in 
the I-V curves in a single-Cgo transistor device, and those 
sub-steps are attributed to the coupling of electron and 
the Ceo-surface vibration mode. In another experiment 
by Leroy et al.^^ the current and the conductance of 
a suspended individual single-wall nanotube device are 
measured, and the phonon-assisted sub-peaks on the two 
sides of the main resonant peak are clearly visible in the 
differential conductance versus gate voltage, which is due 
to the radial breathing phonon mode. On the theoretical 
side, the influence of e-ph interaction on the mesoscopic 
transport is also studied by several groupsj ^^i^^i^^ Many 
interesting results, e.g. the phonon-assisted sub-peaks, 
etc, are first theoretically predicted, and then experimen- 
tally observed. 

In this paper, we investigate the transport behavior 
of a molecular QD system having an e-ph interaction by 
using the scattering matrix method. We focus mainly 
on the transmission phase of the phonon-assisted tunnel- 
ing sub-peaks, as well as the phonon-induced dephasing 
process. The results exhibit that the transmission phase 
6 drops between two adjacent (sub)-peaks and 9 rises 
again near the position of sub-peaks. In particular, the 
characteristic of phonon-assisted tunneling process in the 
transmission phase is much more pronounced and visible 
than these sub-peaks in the conductance. Afterwards we 
discuss the dephasing ratio. While at zero temperature 
and at low bias V (V < loq with loq being the phonon 
frequency), the electronic transport through the molecu- 
lar QD is completely coherent because the electron can 
not absorb or emit phonons under this condition. How- 
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ever, if at non-zero temperature or at a high bias Vbias 
{Vbias > ^o)i the dephasing process occurs. In the Umit 
of high bias Vbias {Vbias ^ (jJo); the dephasing ratio goes 
as square of the e-ph interaction strength A in the weak 
interaction region, but it is hnearly dependent on A in the 
strong interaction region. In addition, we also consider 
an open AB ring device with a molecular QD embedded 
in one of its arms, and find that it is feasible to exper- 
imentally measure the influence of the e-ph interaction 
through the transmission phase. 

The rest of this paper is organized as follows. We in- 
troduce the model and derive the formula of transmission 
amplitude in Sec. II. In Sec. Ill, we present the numer- 
ical results and their discussions. In Sec. IV, we study 
the phase measurement by using an open AB ring device. 
Finally, a brief summary is presented in Sec. IV. Some 
detailed derivation of the transmission amplitude is given 
in Appendix. 



II. MODEL AND FORMULATION 

The system under consideration is a molecular QD cou- 
pled to left and right leads in the presence of a local 
phonon mode, and it can be described by the following 
Hamiltonian: 



H — Ho + Hi, 



where 



Ho 
Hi 



Qr,fc 



Xib^ + b)dU + J2[to.kcl^d + H.c]. 



(1) 

(2) 
(3) 



Here cl^^.{cak) and S{d) are the electron creation (an- 
nihilation) operators in the lead a = L, R and the QD, 
respectively. 6^(6) is the phonon creation (annihilation) 
operator in the QD. Due to large level spacing of the 
molecular QD, only one relevant quantum level Eq is con- 
sidered. The electron in the QD is coupled to a single 
phonon mode wq, and A and tak describe the strength 
of the e-ph interaction and the coupling between the QD 
and the leads, respectively. 

In the following, we apply the S-matrix scattering for- 
malism to derive the transmission amplitude, the trans- 
mission phase and the current. From Hamiltonian (1), 
the S matrix can be written asi ^^'^"^ 

/oo 
dtie'"°''Hie^'"°''e''^^''^ 
-OO 

dtidt2e'""^^HiGrit2-h) 



oo roQ 



— oo — oc 



xHie 



-iHoti -v{\ti\ + \t2\) 



0^ 



(4) 



where the single-particle Green's function operator (t) 
is, Gr{t) = —iO{t)e~^^* . By using the S matrix, the 



final state |/ > can be obtained from the initial state 
|i >, with 1/ >= S" X |i >. Considering an initial state 
\i >= \ei,n,L >, which denotes an electron with energy 
£i in the left lead and n phonons in the QD, the final 
state 1/ > can be expressed as: 

1/ > = 5 X |i >= 5 X \e^,n,L > 

= ^ [r,n{ef,ei)\ef,m,L>+Un{ef,ei)\ef,m,R>] 

m— — OO 

+ 00 

= ^ [rrni£i)S{si+nLL>o- £f -mLL>o)\sf,m,L > 

m— — oo 

+tm.i£i)S{e^ + nujo - £/ - mwo)|e/,m, R >]. 



(5) 



where tm(e) and rm{£) are the transmission amplitude 
and the reflection amplitude with accompanying absorp- 
tion or emission of |m — n| phonons. At zero tempera- 
ture, the phonon number n in the initial state \i > must 
be zero and then tm{£i) can be written as (the detailed 
derivation is shown in the Appendix): 



/ 



tm{£i) = / dej t„i{ef,£i) 



'm 



1=0 



ll{m-l)\ 



j G^iSi - nojo - li^o) 

77,' 



(6) 



Obviously, to{e) describes the amplitude of an elastic tun- 
nehng process which is coherent. While tm(e) {m ^ 0) 
is the amplitude of an inelastic tunneling for emitting 
m phonons. Due to emission of phonons, thus, leaving 
a trace in the QD for the inelastic tunneling process, 
an inelastically tunnelled electron loses its phase coher- 
ence. So at zero temperature the transmission phase shift 
through the QD isj^i^ 



^arg{to{Q)}. 



(7) 



From tm{e), the total transmission probability (including 
the coherent and the non-coherent parts) through the QD 
is Ttotis) = J2m=o ^nd the transmission proba- 

bility of the non-coherent part is Td{e) = J2m=i 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

In this section, we numerically study the transmission 
phase 9 and the dephasing ratio Td/Ttot- In our numer- 
ical calculations, the phonon frequency ljo is set as the 
energy unity (wq — !)• Notice that the main result of 
this paper, Eq. (6), is obtained at zero temperature, and 
thus so are the numerical results and their discussions. 
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However, results at low temperature should be similar 
to that at zero temperature. Fig. 1 shows the coherent 
transmission probability |top [namely, |to(0)p] and the 
phase 9, function of renormalized level Eq, which 

can be tuned by the gate voltage in an experiment. No- 
tice that 1^0 P is proportional to the linear conductance 
G through the QD, G = {e'^ / h)\to\^ . Due to the e-ph 
interaction, several interesting features arc manifested. 
In addition to the main peak related to the single level, 
new satellite sub-peaks appear in the curve of Itop-S'o at 
— Eo = nujQ [n = 1,2,...). The sub-peaks only exist on 
the right hand side of the main peak, and their heights 
increase with the increase of the e-ph coupling strength 
A. The sub-peak at eo = corresponds to the phonon- 
assistcd tunneling process as shown in the inset of Fig. 
1(a), in which an incident electron from the left lead first 
emits a phonon and tunnels to the level Eq, and subse- 
quently reabsorbs a phonon and tunnels forward to the 
right lead. Since this process does not leave a trace in 
the QD, it maintains the phase coherence. Meanwhile at 
zero temperature, there is no phonon in the QD in the 
initial state and the absorption process can not occur, so 
that the satellite sub-peaks only exist on the negative 
side [see Fig. 1(a)]. 

Next, we study the transmission phase 6 that exhibits 
a non-monotonic behavior. Across the main resonant 
peak, continuously rises by a value of tt. This result 
is consistent with the previous theoretical and experi- 
mental findings.— i^i^ Because of the c-ph interaction, 6 
drops between the main peak and the next sub-peak or 
between two adjacent sub-peaks, and rises again across 
a sub-peak, such that a dip appears around — £o = nuj^. 
These dips are much more pronounced than those sub- 
peaks in the transmission probability |to|^. For example, 
for A = 0.8 the 2-nd phonon-assisted sub-peak is so small 
that it is hardly visible [see Fig. 1(a)], however even the 
4-th dip can be clearly seen [see Fig. 1(b)]. The sensi- 
tivity of the transmission phase to the e-ph interaction 
provides a new way to detect the strength of the e-ph 
interaction. 

In the calculations above, the tunneling coupling 
strength F (F = F^ = F^) between the leads and the 
QD is set to be quite weak, F = 0.1 <C wq. With an in- 
crease of F, the phonon-assisted sub-peaks and the main 
peak in the curve of [tol^-S'o gradually merge together and 
become indistinguishable, and the dips in the curve of 9- 
eo are also gradually getting smaller (see Fig. 2). While 
F w Wo (e.g. F = 0.7), all sub-peaks and all dips are 
almost invisible. Consequently, in order to experimen- 
tally detect the phonon-induced dips of the phase 9 or 
the phonon-assisted sub-peaks, the coupling strength F 
should be tuned to be less than u)q/2. In fact, the condi- 
tion F < ti;o/2 is normally satisfied in the experiments 

Let us study the amplitude t„i{e) (m = 1, 2, ...) of the 
inelastic tunneling process. In this inelastic tunneling 
process, an incident electron emits m phonons while tun- 
neling through the QD. However it is prohibited when 
Vbias < rnojQ. At small bias case Vbias < and at zero 



temperature, all inelastic tunneling processes are pro- 
hibited and the tunneling through the QD is coherent. 
On the other hand, with Vbias > i-^q, inelastic tunnel- 
ing processes occur and the tunneling through the QD is 
partly non-coherent. In the limit of larger bias voltage, 
Vbias 2> too, the total dephasing transmission probability 

oo 

Td{e) is: Td{e) ~ J2 Km|^- Fig. 3 shows the dephasing 

m— 1 

transmission probability (or Td/Ttot) versus renormal- 
ized level Eo and the e-ph interaction strength A. While 
without the e-ph interaction (namely A = 0), no inelastic 
tunneling process happens and thus Td = 0. When A 7^ 0, 
the inelastic tunneling process occurs and Td is no longer 
zero. A serial of peaks are exhibited in the curve of Td-So 
and the interval between the two adjacent peaks is loq [see 
Fig. 3(a)]. As is seen from Fig. 3(a), a higher peak must 
correspond to a larger value of A (including the peak at 
£0 = 0). This means that the dephasing probability Td 
monotonously increases with the A regardless of the po- 
sition of the renormalized level Eq. Next, in Fig. 3(b) 
we show the relative dephasing transmission probability 
Td/Ttot versus the e-ph coupling strength A in the reso- 
nant tunneling region (namely Sq = 0). When A is small 
(A < 0.2a;o), the relative dephasing transmission proba- 
bility Td/Ttot increases parabolically with A, but the de- 
phasing probability Td/Ttot is found to increase linearly 
with increasing A between the range of 1 > A > 0.4. For 
large A case (A > 1), Td/Ttot > 0.6 and the dephasing in- 
elastic tunneling processes dominate. In an experiment 
the parameter g = (X/luq)'^ is generally in the range from 
0.1 to 1,— though some special devices— show a big vari- 
able range of g. In this A region, the degree of dephasing 
is linearly dependent on the e-ph coupling strength. 



IV. THE AB RING DEVICE 

In Sections II and III, we only consider a simple device 
consisting of a QD coupled to two leads. However, in a 
real experiment to measure the transmission phase 9, the 
device is an open AB ring with a QD embedded in one 
of the arms.— Therefore, it is of experimental relevance 
to study the open AB ring device in this section. A QD 
is embedded in one arm of the ring, and the other arm is 
the reference arm with the transmission amplitude tref- 
Due to openness of the open AB ring device, the pro- 
cesses of multi-time circling around the ring is negligible. 
Note that only elastic tunneling process to{ei) is in in- 
terference with the reference arm. TAB{£f,£i), defined 
as the probability that an electron of energy Si incident 
from the left lead will be transmitted with energy s / into 
the right lead, can therefore be written as: 

00 

TAB{£f,£i) = ^ S{ei ~£f- mujQ)\tm{£i) + <5m,o e^'^i^e/l^, 

m=0 
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where (jj is the magnetic flux inside the ring. In the ab- 
sence of the reference arm (namely, t^ef = 0), T/iB(e/, e^) 
is reduced toT(e/,£i) = J2Z=o^i'^-i~'^f~'"^^o)\tnii£i)\'^, 
and this result is the same as that in the work by 
Wingreen et alJ^ Using the transmission probability 
TAsi^f, £i), the current flowing through the AB ring issiS 

Jab = ^ J dsi J defTABi£f,£i)fLi£i)[i - fnisf)] 

J de^ j defTAB{ef,ei)fR{ei)[l~ fL{ef)], 

(9) 

where fhisi) = f{e, -iiL) and = f{ef-HB.) with 

the chemical potential fJ-L{B.) ~ ±eVf„;as/2 and /(e) is the 
Fermi distribution function. Finnaly, the differential con- 
ductance Gab can obtained from Gab ~ dJAB/dVbias- 

Based on Eq.(lO), we show the numerical results of 
the current Jab flowing through the open AB ring in 
Fig. 4. Fig. 4(a) and Fig. 4(b) correspond to small and 
large bias voltage cases, respectively. In both cases, the 
phonon-assisted sub-peaks can be seen on the right hand 
side of the main peak, and the sub-peak height increases 
with increasing e-ph coupling strength A. These results 
of current are similar to those in the previous paper.— 

Next, we focus on the differential conductance Gab 
{Gab — dJAB /dVbias ) and its dependence on the renor- 
malized level Sq or on the magnetic flux (f>. In fact, Sq 
and (j) can be well controlled and are continuously tun- 
able in an experiment. The differential conductance Gab 
is always a periodic function of the magnetic flux (f) with 
a period of 27r. Figs. 5(a) and 5(b) show the linear 
conductance Gab at zero bias voltage. Here a series of 
phonon-assisted sub-peaks exhibits in the curve of Gab 
versus eg, similar to that in the transmission probability 
l^oP [see Fig. 1(a)] because of small value of tre/- Be- 
sides, the phonon-assisted tunneling processes can also be 
observed from the amplitude of the Gab oscillation ver- 
sus the magnetic flux (p [see Fig. 5(b)]. While Eq = 
or —1 (i.e. at the main peak or the 1-st sub-peak), 
the AB oscillation amplitudes are quite large since the 
phonon-assisted elastic tunneling processes play a role 
here. But at the position between two adjacent peaks 
(e.g. eq = —0.5 or —1.5), the AB oscillation amplitude 
is quite weak. When a small bias voltage is applied be- 
tween the left and right leads, all peaks in the curve of 
Gab-£o, including the main peak and phonon-assisted 
sub-peaks, split into two [see Fig. 5(c)], and their po- 
sitions are at eq = mujo ± Vf,ms/2. The reason is that 
at these values of Sq ~ mujQ ± Vbias/'^, the renormal- 
ized level Eq is in line with the left or the right chemical 
potential fiL,R = ±Vf,ias/2, or the distance between sq 
and hl,r is just mujQ. The behavior of the conductance 
Gab versus the magnetic flux cf) for small bias is similar 
with that of the linear conductance [see Figs. 5(b) and 
5(d)]. Note that at zero bias or at small bias {Vtias < ^o) 
all tunnclings through the QD are completely coherent, 
and the small amplitude oscillation in Gab is due to the 



small transmission probability |io|^- Finally, we investi- 
gate the large bias case (Vbias > 'j-'o)- At large bias Vbias 
the peaks in the curve of Gab versus Eq clearly split 
into two with an interval of Vbias- Moreover, some ex- 
tra sub-peaks emerge even on the left of the main peak. 
For example, the sub-peak marked by "B" in Fig. 5(e) 
stands on the left of the main peak at Sq = — Vbias I and 
their interval is wq- In fact, this peak is from the inelastic 
tunneling process t\(£) and a phonon is left in the QD 
with an electron tunneling through the dot. Fig. 5(f) 
shows the conductance Gab versus the magnetic flux (/> 
while the level Eq is fixed on the peak positions of Fig. 
5(e). The amplitude of the AB oscillation of the peak 
"B" is very weak, but the amplitudes are quite large for 
other three peaks. This gives a proof that the peak "B" is 
indeed from the inelastic tunneling process and the cor- 
responding tunneling electron losses its phase coherent. 

Since the differential conductance G ab is a periodic 
function of the magnetic flux 0, one can make the Fourier 
expansion: GAB^ffy = G\g + G'j^gCOs{(f) -|- ^o)- Here 
the initial phase 9q is a direct experimentally measurable 
quantity. Let us compare the measured phase with 
the transmission phase 9. Fig. 6 shows the phase ver- 
sus eo for different bias Vbias ■ While at zero bias voltage 
the characteristics of the phase including the phonon- 
induced dips, are completely the same as the transmis- 
sion phase 9 (see Fig. 1). When a small bias voltage 
(e.g. Vbias — 0.2 < Wo) is applied between the two leads, 
the phase 6*0 changes slightly but can still reflect the 
transmission phase 9 quantitatively. So at zero or small 
bias the transmission phase 0, including the intriguing 
characteristics due to the e-ph interaction, can be di- 
rectly observed through the measurement of differential 
conductance versus the intra-dot level. For a molecu- 
lar QD device, the phonon frequency is usually from 
5mev to 35mev) ^^'^^ so the condition Vuas < is easily 
reachable. On the other hand, at large bias case (e.g. 
Vbias = I-Swq), the phase deviates severally from the 
transmission phase 9. 



V. CONCLUSIONS 

In summary, we study the influence of the electron- 
phonon (e-ph) interaction on the transmission phase and 
the dephasing while electron tunneling through a molecu- 
lar quantum dot. It is found that the transmission phase 
versus the intra-dot level exhibits a non-monotonic be- 
havior, and a pronounced dip emerges when the renor- 
malized level locates at the position of the phonon- 
assisted sub-peaks. In particular, phonon-induced dips in 
the transmission phase are much more apparent than the 
phonon-assisted sub-peaks in the conductance. Besides, 
phonon-induce dephasing increases monotonically with 
the e-ph interaction strength A. The dephasing probabil- 
ity Td is proportional to A^ at small A, but T^^ cx A for 
large A. In addition, the open AB ring device is investi- 
gated. At zero bias or small bias, the measurement phase 
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from the differential conductance versus the magnetic 
flux is found to have the same characteristics with the 
transmission phase, including the phonon-induced dips. 



where \tLiR){e)\^ = j:k\tLHB.k)\^Sie - SLk(Rk))-^^ At 
zero temperature, the above equation can be rewritten 
as: 
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APPENDIX 

In this appendix, we present a detailed derivation 
for the transmission amplitude Eq. (6). Be- 

cause those states are normalized according to {s,n,a \ 
e', n', a') = ^n,n'^a,a'<5(e — only the last term in the 
S-matrix [Eq. (4)] contributes to the scattering matrix 
element tm{sf,ei).— Therefore the scattering matrix cl- 
ement t,„(£/,£i) is reduced into: 



ical transformation withii^ H = e/He"'' and s = 
(A/ti;o)(6t — b)d''d,. Under this canonical transformation 
, Hamiltonian (1) becomes: 



H = H, 



H 



ph 



(13) 



where 



a.k 



(15) 



tmisf, Ei) = (ef,m,R \ S \ Si, 0, L) 



-0{t2-h) 



CO />oo 



oo ^ — oo 



where eq = Eo — gtuo is the renormalized level of the 
QD and tak = takX, with g = (A/wq)^ and X = 
exp [— (A/a;o)(5t + b)]. Next we employ the same approx- 
y^gi{e]+mu]o)t2-ieiti (^^^ ^ ^ | /f-|^g-«^f(*2-ti)^^ | g. Q ^ymation as one in the paper (— ), tak ~ tak- Under this 

approximation, e-ph interaction can be decoupled and 
^m{£f,£i) in Eq. (13) varies into: 



To take the change of variables: ti = ti and t = t2 — ti, 
the integration over ti now yields a 5 function of energies 
as 1] 0"'', and trn{sf,Si) changes into: 

(£/,£») = {ef,'m,R \ S \ e^,0,L) 



tm{£f,£i) = -2TTd{t)S{ei- Ef - mujo) 



tBi£f)tl{e,) 



dh 



Aci t 



Tr{b"'{t)d{t)X{t)d^X^ 



— 27r(5(£,; — £f ~ muJo)0{t) 



dti 



Ssit 



-i2i:5{ei — £/ — mwo) 



-f,m,R\ Hie-'"* Hi \ e^,0,L) 



= —2TTS{ei — Ef ~ muJo)0{t) 



dti 



-isit 



EE 

fc' k 



X / dte'-''*G''(t)e-*™(*=-*i\ 



xtRk'tlk{ef,m,R I c^fe,de '"*d''cLk \ ej,0,L) 



= —2TTS{ei — Ef — muJo)0{t) 



dti 



xtR{ef)tl{e^){m \ de~'"'d^ \ 0) 



-2'K5{£i — Ef — mujQ)6{t) 



dte"- 



xtR{ef)tl{e,){0 I ^de-'"'d^ \ 0) 



— 27r(5(ei — £/ — mujQ)9{t) 



dti 



Acit 



xtR{ef)tl{e,){Q\^d{t)d^ 



(16) 



where G"^(0 = -ie(t)Tr^i{d{t)d^ = -ie{t){0 \ d{t)d^ \ 
0) and e-*"(*^^*i) = Trph{h"'{t)X{t)X^. Using 
the method of Feynman disentangling of operatorsji^ 
g-*m(t2-ti) gg^jj obtained: 



-*™(t2-^ti) ^ rrp,,{fe"(t)X(i)Xt} 

= (0 I b"\t)X{t)X^ I 

= e-^"(0 I 6"(i)e^*"*e-''" | Q)ph 

= e-^"e-*™"°*(0 I &"e''*"* | 0)p,, 



(11) 



ml 



(17) 



6 



where m = A[l - e-^'^"(*2-ti)]^ Substituting Eq. (18) into 
Eq. (17), we have: 



tR{£f)tl{£:) 



X / dte^'^'Cit) 



e e 



Te 



-A- 



E 



(-1)"-'™! 



(=0 

°° \2n+m 
: — ^ A -^T 



E 

n=0 



l\{m-l)\ 
I G' (£j - nwo - Zwo), (18) 



where G {E) is the Fourier transform of G {t). Here 
we have assumed the symmetric couphng [ ti^{ei) = 



tR{£i)] and considered the wide-band Umits case, so 
r = 27rti(fl)(e/)t2|.j^^(ei) is independent of the energy 
Ei and £/. In the wide-band hmit the Green's function 
G"'(e) is easily calculated following the standard proce- 
dure^^>2i 



G"'(e) 



£ — £o + iT 



(19) 



From tmi£f,£i), the transmission amplitude tm(ei) can 
be obtained: tm{£i) = J de ftm{£ f T£i)^ and the result is 
given in Eq. (6) in the text. 
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FIG. 1: (Color online) The transmission probability |toP (a) 
and the transmission phase 6 (b) vs. the renormalized level eo 
for the different e-ph interaction strength A with F = 0.1. The 
inset in (a) is the schematic diagram for the phonon-assisted 
tunneling process. 



FIG. 3: (Color online) (a) e-ph coupling strength A depen- 
dence of the dephasing probability, Td- (b) The relative de- 
phasing probability Td/Ttot vs. A. The parameter F = 0.1 in 
(a) and (b). The dot line in (b) is guide to the eye. 



FIG. 2: (Color online) The transmission probability |toP (a) 
and the transmission phase 9 (b) vs. the renormalized level eo 
for the different F with the e-ph interaction strength A = 0.7. 



FIG. 4: (Color online) The current Jab vs. the renormalized 
level £o for the different e-ph interaction strength A with F = 
0.1, <^ = 0, and Uef = 0.1. 
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FIG. 5: (Color online) (a, c, and e) are the differential con- 
ductance Gab vs. the level eo for the bias Vbias =0 (a), 0.2 
(c), and 1.5 (e) at = 0. (b, d, and f) are the differential con- 
ductance Gab vs. the magnetic flux 4> for the bias Vbiaa ~ 
(b), 0.2 (d), and 1.5 (f). The other parameters are F = 0.1 
and tref = 0.1. 



FIG. 6: (Color online) The phase 9o vs. the level eo for 
the different bias Vbias- The other parameters are A — 0.6, 
F = 0.1, and Uef = 0.1. 
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